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Abstract: A review of theoretical models of diffractive structure functions in deep 
inelastic scattering (DIS) is presented with a view to highlighting distinctive fea- 
tures, that may be distinguished experimentally. In particular, predictions for the 
behaviour of the diffractive structure functions F®, ', ^7 i 2 D ( c ^ iarm ) are presented. 
| The measurement of these functions at both small and high values of the variable 

(3 and their evolution with Q 2 is expected to reveal crucial information concerning 
the underlying dynamics. 
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1 Models of hard diffractive structure functions in DIS 



D , It is natural to start with a definition of what we mean by the terms 'hard' and 'diffractive' 



when applied to scattering of electrons and protons. High energy scattering processes may be 
conveniently classified by the typical scales involved. By hard scattering we mean that there 
^ ■ is a least one short distance, high momentum, scale (e.g. high pr-jet, boson virtuality, quark 
mass) in the problem that gives one the possibility of using factorization theorems and applying 
perturbative QCD. In case of diffractive DIS this is the photon virtuality, Q 2 , however this hard 
scale is not necessarily enough and indeed QCD factorization may not even be applicable to all 
hard diffractive scattering in DIS (see |TJ, ^, [3| for discussions and refs). It has been shown to 
be applicable to diffractive production of vector mesons initiated by a longitudinally polarized 
photon [[|. For the time being we will use the definition, due to Bjorken, that a diffractive 
event contains a non-exponentially suppressed rapidity gap. Rapidity is the usual experimental 
variable related to the trajectory of an outgoing particle relative to the interaction point: 
given approximately by rj ~ — ln(tan(0/2)) (in a cylindrical system of co-ordinates centered 
on the interaction point, with the z-axis along the beam pipe and polar angle 6). This rather 
obscure sounding definition results from the fact that within perturbative QCD large rapidity 
gaps (LRG) are suppressed because a coloured particle undergoing a violent collision will emit 
radiation that would fill up the gap. The suppression factor increases with the interval of 
rapidity but it's absolute magnitude for diffractive processes in DIS is uncertain. An additional 
source of rapidity gap supression comes from an overall damping factor associated with multiple 
interactions. The amount of damping is found to be much smaller in DIS than that typical for 
soft processes (e.g. proton proton collisions see ||) making LRG events more likely. 
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Theoretically, for ' diffract ive' electron proton scattering in DIS one must observe the proton 
in the final state. In practice this is very difficult for HERA kinematics since the highly 
energetic scattered proton disappears down the beam pipe in most events. This means that 
the current measurements also contain contributions from interactions in which the scattered 
proton dissociates into higher mass states. This uncertainty is considerably alleviated by the 
advent of the Leading Proton Spectrometer (LPS) which will provide crucial information about 
diffraction (for the first data from the LPS see M). The significance of the difference between 
the experimental working definition of diffraction and the theoretical one is an interesting but 
as yet unresolved problem (it is certainly possible to produce large gaps in rapidity in 'non- 
diffractive' processes, e.g. via secondary trajectory exchanges). 

Such LRG events occur naturally in processes known to be governed by soft processes (e.g 
proton anti-proton scattering at high energies). This is explained naturally in the context of 
Regge theory : at high enough energies one reaches the so-called Regge limit (s>t and s ^> all 
external masses) and all hadronic total cross sections are expected to be mediated by Pomeron 
exchange and to exhibit the same energy behaviour. This expectation is born out by the data 
(see e.g. 0), which shows that a wide variety of high energy total elastic cross sections have 
the same energy dependence which is attributed to the trajectory of the soft pomeron. The 
energy dependence for diffraction in these processes is discussed in e.g. [§]. 

Scattering of virtual photons and protons at small enough x corresponds to the Regge limit 
of this subprocess (s ^ t, s ^> Q 2 , M 2 ioton ). It is natural to ask if the diffractive events observed 
in the DIS sample also exhibit the universal behaviour even though we are now considering off- 
shell scattering for which, strictly speaking, Regge theory does not necessarily have to apply. 
One of the reasons why hard diffraction at HERA at small x is so interesting is that as x 
decreases, for fixed large Q 2 there should be a transition between the hard short distance 
physics associated with moderate values of x and the physics of the soft pomeron which is 
widely believed to dominate at very small x. It is a theoretical and experimental challenge to 
establish whether LRG events in DIS in the HERA range are governed by hard or soft processes 
or whether they are actually a mixture of both. The purpose of this report is to discuss the 
current theoretical models for diffractive structure functions in an attempt to address this 
problem, and, in particular, to outline the benchmark characteristics of the various approaches 
to facilitate the search for appropriate experimental tests. 

In analogy to the total DIS cross section, the diffractive cross section in DIS can be written, 
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2[l + R D (x,Q 2 } x P} t)] 

where D denotes diffraction, R D = F^/(F^ > — F£) and y = Q 2 /sx ; t = is usually assumed 
since the cross section is strongly peaked here. 

Ingleman and Schlein || suggested on the basis of expectations from Regge theory that the 
diffractive structure functions could be factorized as follows: 



F 2 D (x,Q 2 ,x P ,t) = f(x P ,t)F?((3,Q 2 ,t), (2) 

where Q 2 is the photon virtuality, x P is the fraction of the proton's momentum carried by the 
diffractive exchange and t is the associated virtuality, (3 = Q 2 /(M 2 + Q 2 ) = x/x P , with M 2 
the mass of the diffractive system. The last relation for j3 in terms of x, the Bjorken variable, 



is a good approximation but only holds for negligible t and proton mass [10|. Due to lack of 
information on the remnant proton both x P and t can only be estimated indirectly or have to 
be integrated out. 



The 1993 HERA data ITT], [L2|, [L3|] confirmed the presence of events with large rapidity gap 
between the proton direction and the nearest significant activity in the main detector, in the 
total DIS cross section at the leading twist level (i.e. this contribution persisted to high values 
of Q 2 ). These events constitute approximately 10 % of the total sample (compared to ~ 40% 



o 



in photo-production). As has been known for many years and as Bjorken has recently pointed 
out |14| the fact the diffractive cross section is present in the total sample as a leading twist 
effect (i.e. it 'scales') at large Q 2 and small x does not necessarily imply that the mechanism 
that creates these events is point-like. For a careful discussion of the kinematics of hard and 
soft diffraction in a variety of different reference frames see [T5 . 



The observed events were also not inconsistent with the Regge factorization of eq.(|2]). Since 
the cross section had the same power-like x P dependence (in /) over a the wide range of ((3, Q 2 ) 
that were measured it was tempting to postulate that a single mechanism or 'exchange' was 
responsible for these events. The presence of the gap tells us that this object is a colour singlet 
and since the centre of mass energy was very high, the exchanged object became known as the 
'Pomeron'. From this observation it is natural to ask, following ||, if the partonic content of 
this 'particle' may be investigated by changing (3 and Q 2 , with f3 interpreted as the momentum 
fraction of the pomeron carried by the struck parton; / in this picture is interpreted as 'the 
flux of pomerons in the proton'. 

This approach has led to a plethora of theoretical papers in which the parton content of 
the Pomeron at some small starting scale, Qq, is treated in various physically motivated ways 
(relying strongly on Regge theory). The DGLAP jTjJ equations of perturbative QCD (to a given 
logarithmic accuracy) are then used to investigate the evolution with Q 2 of this parton content. 
Formally the use of the DGLAP equations is inapplicable for the description of diffraction 
because the presence of the gap makes it impossible to sum over all possible final hadronic 
states. Their use in this context is at the level of a plausible assumption. In some papers an 
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take the view that that the Pomeron may be 
more like the photon and so can have, in addition, a direct coupling to quarks within the virtual 
photon. Although it is no longer clear once a direct coupling has been introduced whether the 
concept of a Pomeron structure function has any meaning. 



Fits |10| , [171 , Plj to the 1993 data on diffraction reveal a partonic structure that is harder 
(more partons at high j3) than the proton and that gluons contain a large fraction of pomeron 
momentum (up to 90 %) with a large fraction of these at high f3. Clearly in a quantitative 
sense such statements will depend on the physical assumptions used to parameterize the input 
distribution. However qualitatively these statements are reasonable. The paper of Gehrmann 



and Stirling [ID| is particularly useful in discussing Pomeron structure function models in that 
it discusses and compares two models: model 1 which has only resolved component and imposes 
a momentum sum rule on the parton content and model 2 which also allows a direct coupling 
of the Pomeron to quarks. This leads to rather different predictions for the Q 2 evolution of 
these two models (see curves labelled 'GS(I), GS(II)' in fig.(|l|)). Model 1 evolves in a way 
familiar to the evolution of the proton structure function in QCD, i.e. as Q 2 increases there 
is a migration of partons from high to low (3. In model 2, as a result of the direct coupling 
of the pomeron to quarks (at l f3 = 1'), the high (3 distribution is supplemented and, provided 



the direct component is large enough, one expects an increase of parton densities with Q 2 over 
the whole j3 range, which is also an expectation of the boson-gluon fusion model of |22[] (see 

fig-©)- 





Figure 1: Distribution of x P F®(f3,Q 2 ,x P ) as a function of (3 and Q 2 , at fixed x P = 0.005, for 
various models. For key assignments - see text. 

The high gluon content of the pomeron that comes out of the LO QCD factorizable pomeron 
models indicate that the pomeron structure .R-factor, R D (/3, Q 2 ,x P ) = F^(j3, Q 2 , x v )/F^(/3, Q 2 , x P ), 
where F® = F® — F®, may be considerably bigger (R D ~ 1) than that for the proton 
(R p ~ 0(a s )). Clearly in order to provide a theoretically consistent prediction for F® a 
NLO QCD calculation is required. Such a calculation has been performed by Golec-Biernat 



and Kwiecinski [lTj who consider a model with resolved partons in the pomeron subject to a 
momentum sum rule. For high /3, R is small in such models but it can reach 0.5 for (3 < 0.1. 
It has a much softer dependence on j3 than F 2 D (see 'GK' in figs.([jfj)). 
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Figure 2: R D as a function off] andQ 2 for x P = 0.005. The pomeron structure function model 
'GK' differs markedly from the two gluon exchange models 'NZ,BW, LND,RS,BP' at high and 
low (5. Of these, those based on the dipole approach to BFKL, C NZ,BP\ produce markedly 
different (3 andQ 2 behaviour to 'RS,BW,LND\ 



This picture of the pomeron structure function has been discussed in detail elsewhere and 
will not be repeated in further detail here. For a lucid account of this picture and of the 1993 
The latest results from HI IE 



data see 21 



24] on the 1994 data (which has a factor of 



10 increase in statistics and covers a broader kinematical range) suggest that single particle 
factorization no longer holds over the full kinematical range and that particularly for small (3 
it breaks down, i.e. / in eq.(pj) become j3 (but not Q 2 ) dependent. A possible explanation of 
this is that sub-leading Regge exchanges play an important role 0, |2T 



The paper by Ellis and Ross |2bj calls into question the validity of these parton model 
approaches using kinematical arguments concerning the virtuality of the struck parton. They 



stress the importance of measuring diffractive events at high /3 and predict a slow power-like 
increase with Q 2 in this region in contrast to the logarithmic decrease that may expect from a 
naive QCD evolution. 

This and others models are, broadly speaking, similar in spirit to the old aligned jet model 
(AJM), which is a kind of parton model approximation to the wavefunction of the photon (see 



||T4| , [T5|| and refs.), and it's QCD improved formulation (see and refs.). Consider virtual- 
photon proton scattering at high energies (small x) in the proton's rest frame. In this frame the 
virtual photon, whose energy, q , is the largest scale, fluctuates into a qq at a large distances, 
l c = l/2M p x = qo/Q 2 , from the proton. As Ioffe |28 observed many years ago these large 



distances are important in determining the DIS structure functions. For the HERA energy 
range this 'coherence length' can be as large as 1000 Fm. In other words, at enough high 
energies we may consider DIS as the interaction of the quark anti-quark pair with the proton. 
The transverse size of the pair on arrival at the proton is b\ ~ \jk\. 

In the configuration in which one of the quarks carries most of the momentum of the photon 
a large transverse distance develops between the fast and the slow quark by the time it arrives 
at the proton. This large system, in which the pair is initially 'aligned' along the direction of the 
original photon, essentially interacts with the proton like a hadron. This aligned configuration 
gives a leading twist contribution to F 2 and _F 2 D , the latter being interpreted as the fraction of 
events where the produced pair is in a colour singlet state. Since the slow quark is almost on 
shell, the AJM is similar to the parton model and there is no leading twist contribution to Fl 
from this configuration. 

In the configuration in which the momentum is shared more equally the quarks can stay 
closer together in transverse space and may interact with the proton perturbatively. These 
configurations contribute at leading twist to F 2 (x,Q 2 ) and Fl(x,Q 2 ). In the former the inte- 
gration over the momentum fraction leads to the logarithm in Q 2 (coming directly from the 
box diagram). For such small configurations colour transparency phenomena are expected: the 
emission of initial and final state radiation is suppressed p? . 



A semi-classical calculation |29|, |3(| in which the proton is treated like a classical background 
field, leads to results very similar to those of the AJM. Working in the proton's rest frame, one 
considers the interaction of different kinematical configurations of the highly energetic partons 
in the virtual photon with the soft colour field of the proton. These interactions induce non- 
abelian eikonal factors in the wavefunctions of the partons which can lead to diffractive final 



states. In |30[ the addition of gluon to the final state is considered. Leading twist diffractive 
processes appear when at least one of the three partons has a small transverse momentum and 
carries a small fraction of the longitudinal momentum of the proton. The other two partons 
may have large transverse momentum, this means they stay close together as they move through 
the proton, acting effectively as one parton. This high kx jet configuration, produces the only 
leading twist contribution to F£ at this order (which is constant) and InQ 2 enhancement of 
F® ■ This signals that F® also has leading twist contribution from the configuration in which 
all the transverse momenta are small. Several qualitative phenomenological predictions come 
out of this picture. One expects the ratio F 2 D j ' F 2 to decrease like InQ 2 and there to be fewer 
high-py jets in F 2 D than in F 2 (they appear only at order a s in the former). Leading twist 
diffraction appears at order a s in F® which will be dominated by jets. 



Buchmuller and Hebecker present a model of diffraction in DIS based on the dominant 
process being boson gluon fusion, with the colour singlet state being formed by soft colour 
interactions (SCI). The main point is that diffractive and non-diffractive events differ only by 



SCI, the kinematics are expected to be similar since one gluon carries most of the momentum 
of the exchanged system. This idea has also been developed in [j3~T| |32|] which provides a Monte 
Carlo simulation of SCI. 



The simplest QCD model for diffractive exchange is a pair of t-channel gluons in a colour 
singlet state. Such an exchange is a common feature of many models p3| , 
P^ ] and leads to a diffractive structure functions which are proportional to the gluon density 
squared. The dynamical content of these models differ in the treatment of QCD corrections 
and choice of gluon density and will be discussed in more detail below. 

It may be possible to distinguish these models from those in which soft colour interactions 
play a role |35|, |3D| by comparing F 2 D (x, Q 2 , x P ) with F 2 (x P , Q 2 ) for fixed Q 2 and intermediate 
(3. For the latter the following scaling relation is predicted: 



34, 35 



F 2 {x, Q , x% 



~—F 2 (x P ,Q 2 ) 
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(3) 



where C is a constant. 



In where diffraction is governed by two gluon exchange one expects this behaviour to 
be multiplied by a factor x~ x , where A > 0.08 and will depend on k\ (see below). In the dipole 
approach to BFKL [[3^, |39|, in which the dipoles couple via two gluon exchange a similar result 
is expected but with a larger power x~ A , A = a P — 1 = 12a s ln(2)/7r possibly softened by 
inverse powers of logarithms in l/x P |38, 4C]. Of course, in this case, the individual energy 
dependences of F 2 and F 2 D is expected to be a lot harder. 



In the perturbative QCD approach advocated by Bartels and WiisthofF p3f the coupling 
of the pomeron to the hadronic final state can be derived without any additional parameters 
except the strong coupling. The following ansatz is used for the unintegrated gluon density: 



H^k^Ql)^-^-^-^, (4) 

with the effective scale-dependent pomeron intercept (which explicitly, albeit mildly, breaks 
the factorization of eq.(|) since it depends of Q 2 ) a F (Q 2 ) = 1.08 + 0.1 ln[ln(Q7lGeV 2 ) + 3] 
for Q 2 > 0.05GeV 2 and 1.08 below this. This gluon density is then fitted to the available data 
on F 2 . Predictions for the diffractive cross section (which is proportional to [ip(x P , k 2 -,, Qq)} 2 
integrated over kj>) with qq and qqg in final state are then presented over a wide range of (3. 
Now the relevant scale in a F is the virtuality k\. 

In the limit (3 — > 1 the longitudinal contribution, which is formally 'higher twist', is finite 
so is expected to dominate over the transverse part which goes like 1-/5. This highlights the 
fact that the concept of 'twist' must be applied very carefully in diffraction - contributions 
which naively appear higher twist may in fact dominate at high Q 2 in certain regimes. With 
an additional gluon in the final state one finds a (1 — (3) 3 behaviour at large (3. For small (3 this 
configuration dominates and the cross section diverges like l/f3. In summary, a characteristic 
(3 spectrum is found that shows that emission of the additional gluon is bound to the small (3 
region whereas the large f3 is dominated by the longitudinal photon. Numerical results, labelled 
'BW, using the ansatz of eq.(f|) for F®, and R D as function of (3 and Q 2 are shown in figs.([l],0). 

The large mass, small (3 or triple Regge regime (s ^> M 2 Q 2 ^> Aq CD ) has also been 
investigated in detail by Bartels and Wiisthoff (see |42|, |43| and refs). Theoretically the 



emergence of a 4 gluon t-channel state which builds up the large diffractive mass is expected. 
Experimentally, this region is hard to investigate since the requirement of a large mass tends 
to close up the rapidity gap making it difficult to distinguish experimentally from the non- 
diffractive background and also because the diffractive final state may not be fully contained 
in the main detector. This situation is improving now that the first data collected with the 
LPS is becoming available |J. For the purpose of this report we will discuss expectations in 
the not-too-small (3 regime. 



Diehl has calculated the contribution of qq in the final state to the diffraction cross 



section in the non-perturbative two gluon exchange model of Nachtmann and Landshoff [g6|. |44 



Numerical predictions for this model (applicable for not-too-small j3), labelled 'LND', are shown 
in figs.([l]j2|j3|). This model predicts a relatively small contribution of charm in diffraction (less 
than 10 % over a wide range of x w , (3, Q 2 ). 

The high (3, small mass regime of diffraction is considered explicitly in who work in 
co-ordinate space of the transverse distance between the quark and the anti-quark. They claim 
that at high enough (3 (> 0.4) only the qq contributes (in agreement with |[33|| ) and that for 
(3 > 0.7 diffractive scattering from the longitudinal photon dominates for which only small 
distances (6r ~ VQ) contribute. The effective scale of the gluon density relevant to diffraction 
is k^/(l — 13) (see also ^5|) which is clearly hard for high (3. This implies that for high (3 (see 



fig-®) R D becomes greater than unity in sharp contrast to the Pomeron structure function 
model of ]T7|]. For the transverse photon distances of 6r ~ lGeV 2 dominate which is used 
to justify the use of perturbative QCD and the use of evolution equations, using GRV input 
distributions, for the diffractive structure functions. 



The series of papers by Genovese, Nikolaev and Zakharov [[f6|, |47], ^8 
for diffraction inspired by the QCD dipole approach |50 



9[ provides a model 



91 51 1 to the generalised BFKL [52 



equation. In jRJ they strongly reject the factorizable pomeron model and instead suggest that 
a two component structure function for the pomeron with valence and sea partons having dif- 
ferent pomeron flux factors. The absolute normalizations of these components of the diffractive 
structure function are substantially the same as evaluated in 1991 [p9| , before the HERA data 
have become available. In recent papers for this Regge factorization breaking model specific 
predictions for Fi |46| and charm j|7| are given (see 'NZ' in figs.(pI0|]|)). 

The curves, labelled 'RS', shown in figs. (|l],|2],|3|) are from a Monte Carlo simulation developed 
by A.Solano and M.Ryskin, for the dissociation of the virtual photon to two and three jets |53| . 
The formulae used are the same as those in the LMRT model but use a GRV [54| gluon 
distribution and a simplified version of the NLO corrections. 



g based on the QCD 



Bialas and Peschanski also present predictions for hard diffraction [ 
dipole picture of the BFKL equation. In this picture they find that most of the diffractive cross 
section comes from the interaction of qq pairs whose transverse size is of the order of the target 
size as seen by the virtual photon. The perturbative QCD prediction is enhanced by the BFKL 
resummation and by the number of dipole configurations in the initial proton state. In the 
factorized picture they find a strong x P dependence modified by log corrections. They expect 
R D to be a strongly varying function of (3 and to go above unity for large (3. The number of 
diffractive events increases with Q 2 over the whole range. At small (3, i.e. large masses, they 
expect a scaling violation to be similar to that seen in F 2 at small x. Predictions of this model 
for F 2 D and Rjj have been presented recently |55| and are shown, labelled 'BP', in figs.(0,|]). 
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Figure 3: Predictions for the charm content in diffraction, as a fraction on the total diffractive 
sample, as a function of (3 and Q 2 for x P = 0.005. The maximum value of ft reflects the charm 
threshold and increases with Q 2 . 
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The ratio of charm events observed in the diffractive structure function is in principle a very 
good test of the hardness of the processes feeding the cc production. Clearly a measurement of 
the 13 and Q 2 spectra for these charm events will provide a lot more information. 

If hard QCD dominates in diffraction, i.e. the transverse momenta of the qq in the loop are 
large, k 2 -. ~ Q 2 , the relative yield of charm in diffraction is determined by the electric charges 
of the quarks and should be about 40 %. In the Pomeron structure function models of jnj the 
charm contribution comes from boson gluon fusion and is indeed large. Model 2 predicts that it 



should also be large at high (3 in comparison to model 1 (compare 'GS(I)' and 'GS(II)' in fig.(Q)). 
Also since diffraction is a higher twist effect one would expect the jT^charm) /^(total charm) 
to decrease quickly CIS db function of Q . 

In the naive AJM, since the quark transverse momenta are small, one would expect a very 
small charm content. Within the QCD-improved AJM this may be expected to increase with 
Q 2 and for sufficiently high Q 2 the charm contribution to diffraction should approach that 
anticipated from hard physics. 

The early paper of Nikolaev and Zakharov [49~|, predicts that the diffractive contribution 



to open charm is around 1 %. In a recent paper [|4"?fl , they present predictions for the charm 
contribution to diffraction and suggest a very steep rise at small x P strongly breaking Regge 
factorization; at x P = 0.005 this leads to a charm content of about 10 % (see 'NZ' in fig.(^)). 

In a numerical study of the influence of the small fcr region in the BFKL equation, in |56| , it 
is shown that the dominant contribution to diffraction comes from the region of small transverse 
quark momenta, even for large Q 2 . This would seem to favour a small charm contribution in 
this model. 

The LMRT approach ]35| is based on the same Feynman graphs for 7 — > qq and 7 — > qqg 



dissociation as §3], 0, [41], |g and [0, g§ @§ [^ |51J . Q However in the LMRT case the 
most realistic MRS (A') gluon distribution (which fits all the present data) were used and the 
main NLO corrections, including an estimate of the K-factor in the 0(—tt 2 ) approximation, 
were taken into account. Thanks to the large anomalous dimension 7 of the gluon structure 
function g(x, k 2 ^) oc (&f.) 7 a ^ s m &U x = x P the infrared divergence is absent from the fc^-integral 
and, even for the transverse part originated by the light quarks, the dominant contribution 
comes mainly from small distances (see also ||37|| ) and doesn't depend too much on the value 
of the infrared cutoff. This short distance dominance is reflected in the large charm content 
of the Monte Carlo [53[] and of |35| (see curves 'RS' and 'LMRT' in fig.(0), respectively). The 



LMRT predictions are normalised using a phenomenological fit to the '93 ZEUS data and show 
a significant threshold behaviour for /3 approaching the kinematical limit. The sharp increase 
for low values of (3 comes from the inclusion of real gluon emission (see p5||), which is not taken 
into account in LND. 

The measurement of the charm contribution in diffraction, which should be available in 
the near future (at least for D* production [|57|), will certainly help our understanding of the 
interplay of hard and soft physics in diffraction. 
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